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O ■ Abstract 

Recently WMAP has measured the cosmological parameters to a much 
greater accuracy. We analyze the implications of this more precise mea- 
■ surement for supersymmetric dark matter and for the direct detection of 

CN . supersymmetry at accelerators. We consider mSUGRA including also the 

\ hyperbolic branch (HB) in the radiative breaking of the electroweak sym- 

^ ■ metry. On the part of the hyperbolic branch where the lightest neutralino 

is dominantly a higgsino rather than being mostly a bino, the relic density 
constraints are satisfied by coannihilation with the next lightest neutralino 



o 

(N 



' and the light chargino. Including this branch the lightest neutralino mass 

■ satisfies m^o < 1200 GeV for tan/3 < 50. Constraints of 6 — > s + 7, of ^f^ — 2, 

. and of — > are also analyzed. It is shown that the neutralino-proton 



Oh! 



cross section in each case will fall within the reach of dark matter exper- 



p 1' iments. Possibility for the direct detection of supersymmetry is discussed 



in the allowed regions of the parameter space consistent with WMAP con- 



D , straints. A brief discussion of the hyperbolic branch and focus point region 



(HB/FP) is also given. . 



1 Introduction 



Recently the Wilkinson Microwave Anisotropy Probe (WMAP) has measured some 
of the cosmological parameters with significantly greater precision[l, 2]. Specifi- 
cally, WMAP gives the matter density of the universe so that VLmh^ = O-ISS'^q^q^qq 
and gives the baryon density so that fi^/i^ = 0.0224±0.0009, where = Pm,h/ Pc 
where pm,b is the matter (baryon) density and pc is the mass density needed to close 
the universe and h is the Hubble parameter in units of lOOkm/s/Mpc . Assuming 
the difference of the two is cold dark matter (CDM) one finds the CDM density 
in the universe according to WMAP is now given by VtcoMh'^ = 0.1126to;oo9. In 
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this paper we analyze the constraint of the WMAP results for supersymmetric 
dark matter. For the analysis we will focus on the mSUGRA model [3] and analyze 
the allowed range of the parameter space consistent with the WMAP relic density 
constraint. The above requires taking account of the full range of the hyperbolic 
branch of radiative breaking of the electroweak symmetry [4]. The mSUGRA model 
is characterized by the parameters rriQ, mi/2, Aq, tan (3 where mg is the universal 
scalar mass, mi/2 is the universal gaugino mass, Aq is the universal trilinear cou- 
pling and tan/3 is the defined by tan/3 =< H2 > / < Hi > where H2 gives mass 
to the up quark and the Hi gives mass to the down quark and the lepton. In the 
analysis we will also consider the b 57 constraint and the 2 constraint. 

tan/3 in the analysis will range up to values of 50 and it is known[4] that for val- 
ues of tan/3 which are large or even moderately large that radiative breaking of 
the electroweak symmetry lies on the hyperbolic branch. To make the discussion 
clearer we review briefly radiative breaking of the electroweak symmetry and dis- 
cuss how the hyperbolic branch arises in such a breaking. One can illustrate this 
phenomenon analytically for the case when the b quark couplings can be neglected. 
In this case one of the constraints of radiative symmetry breaking determines the 
Higgs mixing parameter /j, so that [4] 



Ciml + Csmf/2 + C'^Al + A/^Lp = + \mI (1) 



Here, 
and. 



' ^1/2 + ^^o|, C'2-C2-\^^ (2) 



1/2 



1 3Do-l 2^ ^ , 

^3 = - t\l C, = -^/, A;.L. = (3) 

A//^ is the loop correction. Ei_2 is as defined in Ref.[4] , t — tan/3 and the 
functions e,f,g,k are as defined in Ref.[5]. Further, Dq — 1 — {rrit/mfY and 
m/ ~ 200sm/3 GeV. 

For small to moderate values of tan /3 the loop corrections are typically small 
and further the renormahzation group analysis shows that C2 > and C3 > 0. For 
such values of tan /3 where the loop corrections have reduced scale dependence one 
finds Ci > independent of any scale choice Q for having the radiative electroweak 
symmetry breaking (EWSB). In this circumstance one finds that the radiative 
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symmetry breaking constraint demands that the allowed set of soft parameters tuq 
and m'l for a given value of fi lie on the surface of an ellipsoid. This condition 

2 

then places an upper bound on sparticle masses for a given value of $ which is the 
fine tuning parameter defined by $ = + ^[4]. This is the ellipsoidal branch of 
radiative breaking of the electroweak symmetry [4]. However, it was found in Ref. [4] 
that for typically larger tan ( ^ 7) when the loop corrections to fi are significant 
along with a significant degree of its variation with the scale Q, the above scenario 
does not necessarily hold. One way to see this phenomenon is to choose a value of 
the running scale Qq at which the loop corrections to fi are minimized. One finds 
then that in some parts of the parameter space where rriQ and mi/2 are relatively 
larger the minimization scale Qq occurs in such a region that it leads to a switch 
in the sign of Ci, i.e. sign(Ci((5o))=— 1- In this circumstance one finds that the 
radiative symmetry breaking condition takes the form 



/2 

rrtT. 

where the sign ± is determined by the condition sign{{^ + \)M]r — C2AI) — ± 
and where 

, _ |($o + i)M|-C-A^| _ |($o + |)M|-C-Ag| 

^ IC3I ^ P \Ci\ 

From the above we see that the presence of the relative minus sign leads to a 
drastically different constraint on the soft parameters due to constraint of the 
radiative breaking of the electroweak symmetry. Here for fixed values of Aq one 

finds that mo and m'l lie on a hyperbola and thus these parameters can get large 

2 

for fixed values of /j, or for fixed values of the fine tuning parameter $. This is 
the high zone of the hyperbolic branch of radiative breaking of the electroweak 
symmetry[4]. Remarkably, the soft parameters can be quite large even while the 
value of $ or /i can be chosen to be significantly small. This is a feature which 
gives a significantly different type of mixing of gauginos and higgsinos than the 
usually explored regions of the minimal supergravity model. In the high zone of the 
hyperbolic branch when mi >> /x, an inversion phenomenon takes place, and the 
neutralino mass becomes essentially /i. The above has a drastic effect on sparticle 
spectrum and on supersymmetry phenomenology which we discuss below. 
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2 Sparticle Spectrum in the Inversion Region of 
the Hyperbohc Branch 



As discussed in Sec.l the constraints on mo and mi for fixed ji for the hyperbohc 
branch are very different than for the usual (eUipsoidal) scenario. Here since mo 
and mi can get large for fixed one finds that the squark and slepton masses get 
very heavy and may lie in the several TeV range (The feature of large mo is shared 
by the focus point region of mSUGRA models[6]). We consider here a specific part 
of the hyperbolic branch where mi >> >> Mz- In this scenario then one finds 
the two lightest neutralino states Xij X2 ^'^'^ fight chargino state xf are essen- 
tially degenerate, each with mass ~ We will call this phenomenon "inversion" 
in that the lightest neutralino switches from being mostly a Bino to being purely 
a higgsino. In fact, this is also the case for the second lowest neutralino and the 
lighter chargino since all of them have a common mass [i to the leading order. The 
degeneracy in lifted when corrections 0(M|/Mi 2) and 0{M'^/fi) are included. 
The remaining sparticle spectrum consisting of quarks, sleptons, gluino and the 
remaining charginos and neutralinos are significantly higher and in principle could 
fie in the several TeV range and perhaps beyond the reach of even the LHC. Thus 
the prospects of observing supersymmetry depends on our ability to observe the 
particles Xi: X2 xf iii addition to the observation of the light Higgs boson. In- 
cluding the lowest order perturbation corrections 0(M|/Mi_2) and 0(M|//x) the 
masses of these three lowest mass states in the inversion region at the tree level 
are given by 



The quantities that are relevant for the observability of these sparticles are the 
mass differences 




(6) 



Thus for > the mass pattern that emerges is 



m^o < rrii < m^o 



(7) 



AM^ = m ± - m o, AM° = 



m^o — m^o 



(8) 
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While 'm^±, m^o, and m^a masses lie in the several hundred GeV to TeV (above 
TeV) range the mass differences AM"*" are much smaller and lie in the range f-10 
GeV. The mass differences can receive loop corrections [7, 8] which can be as much 
as 25% or more. However, these corrections do not modify the general picture 
of this scenario. The above leads to some important constraints on what may be 
observed experimentally. 

3 Coannihilation, relic density, and detection rates 
with WMAP Constraints 

We discuss now the WMAP constraints on SUSY dark matter and also investigate 
if such dark matter will be accessible to direct detection. This issue is of great 
importance as there are on going dark matter experiments [9, 10, 11, 12] and also 
experiments planned for the future[13, 14] to detect dark matter. In the analysis 
we will use a 2a constraint on the WMAP[1, 2] result for CDM, i.e., we take 

n^h' = o.iuet'Zl (9) 

Many interesting theoretical investigations in the analysis of supersymmetric dark 
matter have been carried out over the years [15, 8, 16, 17]. These include inves- 
tigations of the effects of the variations of uncertainties in the relic density and 
wimp velocity on the detection rates[18], effects of nonuniversalities in the Higgs 
sector[19, 20] and in the gaugino sector[21, 22], effects of CP phases[23], and the 
effects of Yukawa unification [24, 25]. More recently the effects of coannihilation 
on supersymmetric dark matter have been analyzed[26, 27, 28, 24, 29, 30, 31, 32]. 
This effect becomes important when the mass of the next to the lightest supersym- 
metric particle (nlsp) is close to the mass of the lightest supersymmetric particle 
(Isp) at the time when the Isp's decouple from the background. In such a situation 
the coannihilation processes involving Isp-nlsp and the nlsp-nlsp annihilation must 
be taken into account. The quantity of interest is the number density ria — J^iT^a 
where a runs over the particle types that enter in coannihilation, and n obeys the 
Boltzmann equation 

— = -3Hn- < aeffv > (n^ - nl) (10) 

where H is the Hubble parameter, uq is the equilibrium number density and (jg// 
is the effective total cross section defined by 

f^eff ^^CrahTan (H) 
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where aab is the annihilation cross section of particle a with particle b, and 
fa = noa/no where noa is the density of particles of species a at equilibrium. 
After the freeze out the nlsp's decay to the Isp and thus n becomes the number 
density of the Isp. It was shown that in mSUGRA one naturally has coanni- 
hilation with the sleptons when the neutralino mass extends to masses beyond 
150-200 GcV with processes of the type xh ^ tl,tZ,th, ^4 ^ and 
^'r^'b. ~^ i'^i'^,'y'-f,'yZ,ZZ,W^W^,hh where / is essentially a r. The above coan- 
nihilation processes extend the allowed neutralino range up to 700 GcV[28]. We 
will show that remarkably the relic density constraints can be satisfied on the hy- 
perbolic branch also by coannihilation. However, on the hyperbolic branch the 
coannihilation is of an entirely different nature. Specifically in the inversion region 
the dominant coannihilation is the Xi ~ xf coannihilation followed by Xi ~ X2 
coannihilation, and by Xi ~ Xi aiid by xf ~ X2 coannihilations. Some of the 
dominant processes that contribute to the above coannihilation processes are [33] 

Xlxt, X2Xt ^ eiUi, AW+, ZW+, W+h 

xtxT, X1X2 ^ UiUi, ddi, W+W- (12) 

Since the mass difference between the states Xi ^"^^ Xi is the smallest the XiXi 
coannihilation dominates. 

In the analysis we include the b ^ constraint [34] and the gfj,—2 constraint [35] . 
The constraint arising from — > for large tan (3 is also discussed. The anal- 

ysis of Ref.[36] gives two estimates for the difference a^^^ — a^^: These are [I] a^^^ — 
af^ = 1.7(14.2) X 10-i"[37, 36] and [II] a^^^ - af^ = 24.1(14.0) x 10-i°[38, 36]. 
These estimates also include corrections from scalar mesons to the muon anomaly 
computed in Ref . [36] . Estimate [I] corresponds to essentially a perfect agreement 
and does not put any effective upper limit constraints on the parameter space. In 
our analysis we consider a l.ba range around the central value of estimate [II], i.e., 
we choose 3.1 x 10"^° < {a^^^ - af^) < 45.1 x 10"^°. We attribute the difference 
to supersymmetry[39]. In Fig. 1(a) we exhibit the allowed parameter space in the 
rriQ — mi plane which satisfies the relic density constraint consistent with Eq.(9) 
for the case tan/3 = 10 and /i > 0. The filled dark circles indicate the regions 
which are consistent with the relic density constraints. We note that this region 
includes a lower branch which is the conventional branch where the relic density 
constraints are satisfied due to coannihilation with staus. For the case of Fig. 1(a) 
this extends to mi of about 800 GeV and m^o of about 300 GeV as can be seen 

2 Al 
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more clearly from Fig. 1(b). However, there is also an upper branch where the 
allowed values of mi consistent with relic density run up to the upper limit chosen 
i.e. 10 TeV. The corresponding neutralino mass, however, runs up only to 1200 
GeV because of the phenomenon of inversion discussed in Sec. 2. As can be seen 
from Fig. 1(b) relic density constraints consistent with the WMAP constraints can 
be satisfied in the inversion region for significantly large values of the neutralino 
mass and values of nio up to 16 TeV. The phenomenon of inversion can be seen 
more clearly in Fig. 1(c) where points consistent with the WMAP constraints are 
exhibited in the mi —m^o plane. The imposition of the — 2 constraint [II] elimi- 
nates all of the inversion region and much of the remaining region of the high zone 
of the hyperbolic branch. However, essentially all of the region allowed by the relic 
density constraints is valid if we consider the (7^ — 2 constraint [I]. In Figs. 2(a), 
2(b), 2(c) we give an analysis similar to that of Figs. 1(a), 1(b), and 1(c) except 
that tan/? = 30. Similarly in Figs. 3(a), 3(b), and 3(c) we give an analysis similar 
to that of Figs. 1(a), 1(b), and 1(c) except tan/5 = 50. For the cases of tan/? = 30 
and tan /5 = 50 the h sj constraint is also displayed. In these cases the region 
below the curves labelled b— > S'j is the disallowed region. 

For large tan/? the constraint from B'^ is also of interest[40, 41]. In 

the standard model the branching ratio for this process is B{Bg —>■ ii^ = (3.1± 
1.4) X 10~^ {Vts = 0.04±0.002) while the current hmit from experiment is B(Bg — > 
IJi^jjr) < 2.6 X 10~^. The current estimates are that RUNII of the Tevatron will 
eventually increase the sensitivity for this process to the limit 10~^[41] which still 
falls short of reaching the branching ratio for this process in the standard model. 
However, it turns out that in supersymmetry this branching ratio is dominated 
by the so called counterterm diagram and the contribution from this diagram 
gives the branching ratio a dependence on tan /3 of tan^ /3 for large tan /?. As a 
consequence the B{B^ fJ^^fJ^ ) branching ratio in supersymmetry can get larger 
than the standard model value by as much as a factor of 10^ which brings it within 
reach of RUNII of the Tevatron. However, the B{Bg Ai^/^~) branching ratio in 
supersymmetry is very sensitive to the sparticle spectrum and falls sharply as the 
sparticle spectrum becomes heavy. In Fig. 4 we give a plot of the -5(5° Ai^A'- ) 
constraint in the mo — mi plane. We find that the current experimental constraint 

_ 2 

on B{B'^ f^^fJ' ) does not eliminate any relevant part of the parameter space 
while -B(-B° Ai^/U") — 10~^ can explore the parameter space in mo up to 700 
GeV and in mi up to about 500 GeV. This mass range is far too small to have 
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any influence on the region of tlie liyperbolic brancli we are focussing on in this 
analysis. For this reason tbis constraint is not very effective in the present analysis. 

A quantity of great interest is the spin independent neutralino-proton cross 
section a^Op{SI) on which experimental limits exist from the current dark matter 
experiments so that a^Op{SI) < 10^'^^ cm^. In Fig. 5(a) we give a plot of cr^Op^SI) 
for tan/? = 10 and /i > 0. In Fig. 5(a) the lower rapidly falling curve that 
terminates at m^o = 300 GeV is the branch on which staus coannihilation occurs. 
The upper curve arises from the low zone of the hyperbolic branch while the patch 
to the right is the one that arises from the inversion region of the hyperbolic 
branch. For values of ncutralino masses below 300 GcV the a^Op{SI) cross section 
arising from the upper curve in Fig. 1(a) is much larger than the one arising from the 
lower branch where the relic density constraints arc satisfied due to ncutralino-stau 
coannihilation. We also note that in Fig. 5(a) the patch to the right indicates that 
the scalar cross sections are quite significant even though one is in the inversion 
region. Thus although the direct detection of supersymmetry in the inversion 
region is more difficult, the neutralino-proton scalar crosss are still substantial. In 
the future dark matter detectors[13] will be able to achieve a sensitivity of up to 
10~^^ cm^. We note that a significant part of the parameter space of Fig. 5(a) will 
be probed by these detectors. In Fig. 5(b) we give a plot of the spin dependent 
neutralino-proton cross section a^Op{SD) for tan/? = 10 > 0. A comparison of 
Fig. 5(a) and Fig. 5(b) shows that the spin dependent cross section is typically 
much larger than the spin dependent cross section by 3-4 orders of magnitude. A 
similar analysis for the case tan/3 = 30 is given in Fig. 6(a) and Fig. 6(b) while 
for the case tan/3 = 50 is given in Fig. 7(a) and Fig. 7(b). The conclusions for 
these cases are very similar to the conclusions drawn from Fig. 5(a) and Fig. 5(b). 
Based on these analyses one finds that for tan /3 < 50, the neutrafino mass range 
consistent with the WMAP constraints on the branch corresponding to neutralino- 
stau coannihilation is m^o < 500 GeV and m^o < 1200 GeV for the high zone 

Xi — Xi — o 

of the hyperbofic branch where the relic density constraints are satisfied due to 
coannihilation with the next to lightest neutralino and the light chargino. These 
constraints remain intact under the imposition of the gn — 2 constraint [I] but the 
constraint arising from the inversion region of the hyperbolic branch is removed 
by imposition of the — 2 constraint [II] . 



8 



4 WMAP Constraints and Discovering SUSY at 
accelerators 



The analysis of Sec. 3 shows that Eq.(9) constraints the parameter space very 
strigently. For the usually explored parameter space of minimal supergravity where 
relic density is satisfied in the region of neutralino-stau coannihilation as well as in 
the low zone of the hyperbolic branch with a moderate amount of higgsino in the 
Isp (i.e. without inversion), one finds that the neutralino mass now has an upper 
limit of about 500 GeV for tan /3 < 50 and mo lies in the few hundred GeV range. 
For this case the corresponding sparticle spectrum should all be accessbile at the 
LHC and perhaps some of it may be accessible at RUNII at the Tevatron. Also 
there are some interesting signals for this branch at the NLC[42]. However, on the 
inversion region of the hyperbolic branch of radiative breaking of the electroweak 
symmetry, mo and mi can get as large as 10 TeV or even higher. In this case 
the squarks and the sleptons would he in the several TeV region and hence they 
would be beyond the reach of even the LHC. The light particles in this scenario 
will be the two lightest neutralinos and the light chargino. However, the signals 
for their detection would be significantly different than for the normal scenarios. 
Specifically, in the inversion region of the hyperbolic branch the mass differences 
among Xi, xt ^^^^ X2 so small that the usual signals discussed for the detection 
of supersymmetry would not apply[43]. 

Situations of the type above have been discussed before in Ref. [44] in the con- 
text of string models and in Ref. [45] in the context of Wino Isp scenarios while 
the experimental search for charginos mass-degenerate with the lightest neutrali- 
nos has been analyzed in Ref. [46]. Here the mass scales are significantly different. 
Thus, for example, in the analyzes of Ref[45] the mass difference of the chargino 
and the nearly degenerate neutrahno is in the range of O(IOO) MeV which al- 
lows for charged particle tracks in the detector of the order of few centimeters 
arising from the decay of the chargino to neutralino such as Xi ~^ Xi^^^h 
xt — Xi^'^Ti"'". In the present scenario the chargino and neutralino masses are 
in the several hundred GeV to 1-2 TeV range and their mass difference lie in the 
range of 1-10 GeV. The mass differences are such that the chargino will always 
decay in the detector and the track length will be too small to be visible. Further, 
the conventional trileptonic signal [47] would yield leptons with energies only in 
the few GeV region to provide a useful signal at the LHC [48]. In Ref. [44] it is 
argued that charginos nearly degenerate with neutralinos may be observable in 
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e~^e~ colliders via observation of hard photons in the process e~^e~ IXiXi ■ 
However, a more detailed analysis for the detection of supersymmetry in collider 
experiments is needed for the scenario discussed here. On the whole, the prospects 
for the detection of SUSY signals at accelerators in this scenario look difficult. On 
the other hand quite interestingly this scenario docs provide a sufficient amount 
of dark matter to populate the universe and a part of the parameter space of this 
branch does yield spin independent neutralino-proton cross sections which lie in 
the range of observability of dark matter detectors. We emphasize that much of 
the high zone of the hyperbolic branch and specifically all of the inversion region 
on the hyperbolic branch can be eliminated if the — 2 constraint [II] holds. 
However, the high zone of the hyperbolic branch would not be significantly con- 
strained if the — 2 constraint [I] holds. This points to the importance of getting 
an unambiguous determination of the leading order (LO) hadronic correction to 
5m -2. 

We comment now briefiy on the relation of the hyperbolic branch to the fo- 
cus point region[6]. As discussed in Sec.l we showed that one can find solutions 
to radiative breaking of the electroweak symmetry where mo and mi/2 can get 
large while /i remains fixed and relatively small. These solutions constitute the 
hyperbolic branch. A part of this region also includes the so called focus point 
region. Thus the focus point region is limited to relatively small values of mi/2 and 
consequently mo is also limited from getting very large because of the radiative 
symmetry breaking constraint relative to the case of the hyperbolic branch. Thus 
the focus point region (FP) is truly a subset of the hyperbolic branch (HB). A fur- 
ther discussion of this point can be found in Ref . [49] where the acronym HB /FP 
is used to describe this region. 

5 Conclusion 

The recent WMAP determination of the cosmological parameters, specifically 
ftmh^ and Jlfo/i^, to a much better accuracy than earlier determinations has im- 
portant consequences for the observation of supersymmetric dark matter and also 
for the direct detection of supersymmetry. In our analysis we have identified the 
difference Q^^^ — Q^/i^ as arising from relic neutralinos and analyzed this possibil- 
ity within mSUGRA. One finds that for the region of the parameter space where 
the relic density constraints are satisfied due to the neutralino-stau coannihilation, 
the neutrafino mass fimit is now reduced to m^o < 500 GeV for tan 3 < 50. The 
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spectrum in this case will all be accessible at the LHC with the possibility of some 
sparticles also being accessible at RUNII of the Tevatron. Also some interest- 
ing signals may arise in this case at the NLC. On the high zone of the hyperbolic 
branch including the inversion region, the WMAP constraints are satisfied remark- 
ably to a very high value of the neutralino mass, i.e., up to m^o < 1200 GeV for 
tan/5 < 50. The satisfaction of the relic density even for such large neutralino 
masses comes about because of coannihilation processes exhibited in Eq.(12). As 
discussed in Sec. 2, in the high zone of the hyperbolic branch itlq and m^o can get 
very large and some of the sparticle spectrum may lie outside the reach of even the 
LHC. Thus the squarks, slcptons and ghiinos may be too massive to be accessible 
even at the LHC. Thus the direct observation of SUSY would be very challenging 
if the inversion region of the hyperbolic branch is realized. In this region the only 
light particles, aside from the light Higgs boson /i", arc the sparticles xfjXijX2- 
The mass splittings among them are typically 0(10) GeV and thus their detection 
poses a challenge. Luckily much of the hyperbolic branch and all of the inversion 
region of the hyperbolic branch can be eliminated by a (7^ — 2 signal. This is what 
happens when we impose the Qfj, — 2 constraint [II]. However, imposition of the 
Q/j, — 2 constraint [I] essentially leaves all of the region of the hyperbolic branch 
including the inversion region intact. This points to the need to achieve an unam- 
biguous gfx — 2 constraint by reducing the errors in the leading order (LO) hadronic 
contributions. We also computed the spin independent neutralino proton cross sec- 
tion a^Op{SI) and found that it lies in the range 10~^^ — 10~^^cm^. A significant 
part of this range will be accessible to the future dark matter experiments[13, 14]. 
Implications of WMAP constraints for supersymmetry have also been reported in 
Ref.[50]. 
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(a) A plot in the mo — mi plane of the allowed region consistent with electroweak 
symmetry breaking and the WMAP relic density constraints for the mSUGRA case. 
The input parameters are Aq = 0, tan f3 = 10, /Lt > and the relic density constraint im- 
posed is of Eq.(9). The white region is the parameter space allowed by the electroweak 
symmetry breaking constraints while the shaded region is disallowed. The filled circles 
denote the region allowed by the relic density constraint. The filled circles just below 



the upper shaded region arise from the hyperbolic branch, 
the black line. 
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(b) A plot in the toq 



plane 



of the allowed region represented by 
black circles consistent with elec- 
troweak symmetry breaking and 
WMAP relic density constraints of 
Eq.(9) for the mSUGRA case in- 
cluding the parameter space on 
the hyperbolic branch. The in- 
put parameters are the same as in 
Fig. 1(a) 



(c) Same as Fig. 1(b) except that the plot is 
in the mi — m^o plane. 
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Figure 1: Relic density constraint and neutralino mass range for tan/? = 10 
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(a) Same as Fig. 1(a) except tan/3 = 30. b ^ s + j contour and a^'^^^ contours are 
also shown. 
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(b) Same as Fig. 1(b) except (c) Same as Fig. 1(c) except tan/3 — 30. 

tan/3 = 30. 

Figure 2: Relic density constraint and neutralino mass range for tan (3 = 30 
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(a) Same as Fig. 1(a) except tan/3 = 50. 5 — > s + 7 contour and a^'^^^ contours are 
also shown. 
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(b) Same as Fig. 1(b) except 
tan/3 = 50. 



(c) Same as Fig. 1(c) except tan/3 — 50. 



Figure 3: Relic density constraint and neutralino mass range for tan (3 = 50 
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(a) A plot of the neutralino-proton spin independent cross section a-^a^p{SI) vs the 
neutrahno mass for the allowed region of the parameter space for all the same input 
parameters and constraints as in Fig. 1(a) 
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(b) A plot of the neutralino-proton spin dependent cross section a^Op{SD) vs the 
neutralino mass for the allowed region of the parameter space for all the same 
parameters and constraints as in Fig. 1(a) 



Figure 5: Spin Independent and Spin Dependent Cross Sections for tan P 
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(a) Same as Fig. 5(a) except tan/3 = 30. 
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(b) Same as Fig. 5(b) except tan/3 = 30. 
Figure 6: Spin Independent and Spin Dependent Cross Sections for tan /3 = 30 
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(a) Same as Fig. 5(a) except tan/3 = 50. 
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(b) Same as Fig. 5(b) except tan/3 = 50. 
Figure 7: Spin Independent and Spin Dependent Cross Sections for tan f3 = 50 
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